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Q^. Abstract 

^ : 

i-^i We study the solutions of the Dirac equation in the adjoint representation(gluinos) 

in the background field of SU(2) unit charge calorons. Our solutions are forced to be 
antiperiodic in thermal time and would occur naturally in a semiclassical approach 
to TV = 1 Super-symmetric Yang-Mills theory at finite temperature. 
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1 Introduction 



In this paper we derive analytic expressions for the finite temperature gluino 
zero modes of the Dirac operator in the background field of the Q = 1 
SU(2) calorons. These are self-dual configuration in R^ x 5*1 including the 
well known Harrington- Shepard (HS) solution [Ij as well as the non-trivial 
holonomy calorons [2]- [I]- The periodicity in one direction, to be referred as 
thermal-time, occurs naturally in a path-integral approach to finite tempera- 
ture Yang-Mills theory and, with the inclusion of spinor fields in the adjoint 
representation (gluinos), in its minimal supersymmetric extension. Calorons 
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are thus the natural objects to be considered in a semiclassical approach to 
these theories at finite temperature. They smoothly interpolate between in- 
stantons and BPS monopoles at zero and high temperature respectively [2j-[5j, 
providing a very interesting link between them. One of the required ingredi- 
ents for such semiclassical analysis is the knowledge of fermionic zero modes 
in the background of the caloron field. Although those in the fundamental rep- 
resentation of the gauge group have been known for quite some time [6]-|10j. 
this is not the case for the gluino zero modes. They have been derived only 
recently by two of the present authors [TT], and just for the case of periodic 
boundary conditions in 5*1. These are the relevant modes for supersymmet- 
ric compactifications but not what is needed when studying Af = I SUSY 
Yang-Mills fields at finite temperature. Antiperiodicity in thermal-time has to 
be required in that case. The goal of this letter is to obtain analytic expres- 
sions for the antiperiodic solutions, derived here for the first time even for the 
trivial holonomy, Harrington-Shepard, case. This requires a different approach 
than the one employed in |TT] which was based on the relation between zero 
modes and self-dual deformations of the gauge field, providing only periodic 
solutions. 

The paper is organized as follows. In section [2] we will describe the strategy 
followed to obtain the antiperiodic zero modes and present the analytic ex- 
pressions for the solutions. In section [3] we analyze their properties in several 
relevant limits, paying particular attention to the one in which the caloron dis- 
sociates into a pair of static BPS constituent monopoles. The trivial-holonomy 
HS zero mode solution and the equal mass constituent monopole cases are also 
discussed in some detail. Conclusions and a brief summary of results are pre- 
sented in section m 



2 Formalism 



As mentioned previously, our goal is that of solving the massless covariant 
Dirac equation in the adjoint representation of the group 

^^ = (1) 



in the background field of a Q=l caloron |2]-[1]- This problem has been par- 
tially addressed in Ref . [TT] . The approach that was followed in that paper was 
based on the well-known relation between self-dual deformations of the gauge 
field and the zero-modes of the Dirac operator in the adjoint representation. 
However, the solutions obtained in this way are periodic in thermal-time with 
the same period P (to be taken equal to 1 in what follows) as the gauge field 
itself. Thus, a different strategy has to be set up to derive the antiperiodic 



modes relevant for finite temperature. In what follows we will present the ba- 
sic idea behind our procedure and the results obtained with it. In all technical 
aspects we will rely strongly in the notation and derivations done in Ref. jllj . 

The observation that leads to our solution is the fact that antiperiodic solu- 
tions turn out to be periodic in the double period. Thus, the method of attack 
developed in Ref. |1^ for periodic zero-modes can be carried over if the whole 
problem is seen as living in this duplicated space-time. This replica trick has 
been used by some of the authors in other works p^flS] and is an important 
source of information when dealing with periodic gauge fields. In our case, the 
problem becomes that of finding self-dual deformations of the Q = 2 caloron 
obtained by the replica procedure. Notice that the topological charge is 2 in 
this case, so we expect 4 (CP-pairs) of self-dual deformations. Since the gauge 
field is periodic in the original period, they can be split into those which are 
periodic and those that are antiperiodic in the original period. The former were 
studied in our previous paper and correspond to the ordinary deformations of 
the Q = 1 caloron. Since there are 2 pairs of those, which are periodic in the 
small torus, we expect to find two pairs of antiperiodic zero-modes. Unfortu- 
nately, although some particular solutions are known [14J, there is no analytic 
general expression for the Q = 2 caloron which would reduce the study of 
deformations to the differentiation of the general solution with respect to the 
parameters of the moduli space. In this paper we will thus follow an alterna- 
tive strategy. Incidentally our results could well prove useful in achieving the 
goal of obtaining the most general Q = 2 caloron solution. 

The general formula relating deformations to zero-modes in the adjoint rep- 
resentation is: 

v[/ = iM^7^(I±75)^, (2) 

with the + or — sign depending on whether the solution is self-dual or antiself- 
dual. V is an arbitrary constant spinor and hence, the zero-modes that we are 
looking for can be arranged into two-dimensional complex vector spaces. These 
spaces are generated by any solution "^ and its euclidean CP transform 

^ — , ^c = 75^^* . (3) 



Our formula can easily be shown to satisfy the Dirac equation provided the 
deformation satisfies the background Lorentz gauge condition. 

D^6A^ = Q. (4) 



Using the general ADHM construction one can obtain formulas for the self- 
dual deformations in terms of those for the Nahm-ADHM data. The ADHM 



construction for SU(2) tells us that a self-dual gauge field can be constructed 

as [TS] 

A,{x) = ^ {u%uy , (5) 



where m is a vector in quaternions, F = 1 + u'^u, and the prime denotes 
the traceless part. The vector u is obtained as the solution of the following 
equation 

(A^ - x^ct^)m = q, (6) 



where the quaternionic matrix A'^ and the vector q are x-independent. We 
have introduced the Weyl matrix W^ = (I, ir) whose adjoints are a^ (r are 
the Pauli matrices). In the proof of self-duality one must demand that the 
following matrix 

R = {A^ - x^a^) {A - Xf,a^) +q®q\ (7) 



is real and invertible. For our purpose it is interesting to write down the ex- 
pression of the adjoint zero-modes in terms of the deformations of the ADHM 
data 

5A^ = — {6q^ - uUA)a>'a''d^uj + h.c. , (8) 



where uj = R~^q. To guarantee that the deformations Sq and SA provide 
a self-dual deformation SA^ satisfying the background field gauge condition, 
one must impose certain conditions. These are best expressed in terms of the 
matrix with quaternionic entries JF = M'^SM = J-'f^a^, where M''' = (g. A'' — 
d'fj,Xfj,). The condition then reduces to the hermiticity of ^^ (JF^ = JF^). 

The previous formulas apply for Q = 1 calorons by extending the vector q 
to become a delta-like functional over the periodic functions in one-variable 
z, while 27riA is a covariant Weyl operator with respect to a 1-dimensional 
abelian gauge field A^{z), the Nahm-dual gauge field. After suitable rotations, 
translations and gauge transformations the caloron Nahm data can be taken 
to be ^ 

q^'\z) = piP+ 5iz - 5,) + P_ Siz + 5^)) , (9) 



where P± = (l±r3)/2. The parameter 6i parametrizes the holonomy, becoming 
trivial for and |. Without loss of generality we will assume in what follows 



that Si < S2 = ^ — Si. In the previous formula the delta functions have to be 
taken as periodic functions in z with unit period. The Nahm-dual gauge field 
of the caloron is given by 

i(°)(z) = -27rS,s{Xlxiiz) + Xlx2{z)) , (10) 



where X3 is the position of the ath constituent monopole on the z-axis. They 
can be obtained from the relations irtiXl + m2X| = 0, and X| — X^ = vrp^, 
where ma = 4:7cSa are proportional to the constituent monopole masses. The 
function Xi is the characteristic function of the interval [— (5i,(5i] and X2 that 
of its complementary. 

Eq. (ITj) implies that the Nahm-dual gauge field is self-dual at all but a finite 
number of points. Eq. Q) is then the solution of the Weyl equation except at 
those isolated points. As we will see later the conditions on the deformations 
SA that enter Eq. ([8]), are precisely equivalent to requiring that SA^^ is again a 
self-dual deformation satisfying the background gauge condition. Thus, they 
can be obtained as the solution of the adjoint Weyl equation of the Nahm-dual 
field, up to delta functions. 

Now we should apply this scheme to the replicated caloron taken as a self-dual 
solution in the double torus with period 2/? (remember j3 is fixed to 1). Since 
this caloron now has charge Q = 2 its corresponding Nahm-dual gauge field 
is now a matrix. Using the general construction of Nahm-dual replicas given 
in [12] we obtain 

iH(,)_('4"'W ^ (11, 



where A^^\z) is the Nahm data of the ordinary caloron, and A^{z) is the 
Nahm data of the replicated caloron. 

One may now wonder which is the corresponding q for such a replica solution. 
We will argue that the solution is actually given by 



^"(^)=(.(ou.Vu)^ (12) 



:^ + 2. 



Notice that each of the components of q and A are periodic with unit period, 
but the whole set is periodic with period 1/2 with a twist matrix given by ti: 

A^{^ + l) = nA^{z)ri. (13) 



The quantity q transforms by periodicity as follows: 

q'^iz^ 1/2) = nq\z). (14) 

From here it is possible to use the general formulas of the ADHM construction 
to verify that indeed we obtain a replicated solution. In particular we have 
that u^{z) is given by: 



Now 



:^R 



dzu'^\z)u'^{z) = [dzu^''^^iz)u^''\z), (16) 



which coincides with F — 1 for the caloron. The replicated gauge potential 
follows from 



^?(^) = ^ J dz{u^Hz)d,u^{z)y = Af{x) . (17) 





One might wonder whether the choice of q and A^ are consistent with the 
condition of self-duality in Nahm-dual space, namely that R = M'^M com- 
mutes with the quaternions. To verify that this is so, one must realize that the 
condition of self-duality should hold only in the domain of these operators. 
These are two-component vectors il){z) satisfying 

^{z + ]^) = ni,{z). (18) 



Thus, they should be of the form 



\z 



Thus, gg^ acting on this vector yields: 



^2 I P+(t>{5i)5{z - 5i) + P-(t>{-5i)5{z + 5i; 



^^ ( P+0(<5i)<5(z + (52) + P_0(-(5i)(5(z-(52^ ' • ^^°^ 



The imaginary part of the upper component coincides with 

p^r,{5{z-5,)-5{z + 5{))(^{z), (21) 

which is what is needed to cancel the self-duahty violation. 

Now we proceed to study the self-dual deformations of this replicated caloron 
satisfying the background field condition. We will make use of our general 
formula Eq. ([H]). The conditions following that equation when translated to 
our case become 



R ,uR} _ 4^2 • L,R,JR _ K„R„]R\ ^ 



D^« ^ r^ _ ,^^[A« V;«] = A7:H g«<5gf - 6q^ql^) cy,a^ , (22) 



where t/''^ = bA^a^ = —6A^/{2'n-). The quantities 6ql} are two-component 
column vector whose elements are linear combinations of delta functions with 
complex coefficients. The holonomy fixes that the argument of the delta func- 
tions must he z±6i and z ±Si + ^. Notice that, as anticipated previously, up 
to the delta functions in the right-hand side, the equation adopts the form of 
the Weyl equation for adjoint zero-modes in Nahm dual space. 

Our next step will then be that of finding the solution of Eq. (1221) . Notice that 
both ip^ and 6q^ = Sq^a^ are the unknowns. Without much effort one can 
demonstrate that given a solution one can obtain other solutions by the opera- 
tion tp^ —^ "ip^Q, Sq^ -^ Q'^Sq^, with Q an arbitrary constant quaternion. This 
transformation is associated to the double degeneracy of adjoint zero-modes. 
We must also point out certain subtleties necessary to understand Eq. fl22l) 
and their solutions. The main idea is that the equation must be understood 
as one relating two operators acting on the space two-component functions 
of the form Eq. ( TT9l) . The right-hand side of Eq. ( l22l) acts by multiplication. 
Thus, the left-hand side must be equivalent, when acting over our space of 
functions, to the multiplication by a linear combination of delta functions. 
This imposes non-trivial conditions on the form of Sq^. In what follows we 
will give the possible values for 6q^ that follow from the previous analysis, as 
well as the resulting form for the equation for ip^, skipping all the details of 
the derivation. 

Before showing the equations, we recall that ip^ is a 2 x 2 matrix in (Nahm- 
dual) colour space 

i.^iz) = ("^^fl ^-;^|). (23) 

The boundary conditions specify that it is enough to know the form of ipu and 
ipi2 (the other components can be obtained by translating in z by 1/2). The 



equations for ipu coincide with those for the Q = 1 caloron, and therefore can 
be associated with deformations that are periodic in time. Thus, our sought 
time-antiperiodic zero-modes should follow from the equation 



d,^i2 + r3(Ai)^i2 = 47rV{ P+ {6{z - 6i) - S{z - 62)) 

+P-{5iz + 52)-5{z + 5,))}Q, (24) 

where the function AA is given by 

Ai ^ Af\z) - if (^ + 1/2) = 27rV(x(-<5i, 6^) - xiS^, 1 - S^)) . (25) 

The arbitrary quaternion Q reflects the degeneracy of solutions mentioned 
earlier. Keeping that in mind one only needs to solve the equation for Q = 
and Q = 1. A particular solution is all that is needed, since the general solu- 
tion can be obtained by linear combinations of these ones with quaternionic 
coefficients. The counting matches the predictions of the index theorem. As 
for the periodic case there are essentially two CP-pairs of zero-modes. 

After these considerations we proceed to show the two particular solutions 
that we will need. The first one corresponds to the inhomogeneous equation 
{Q = 1) and is given by 

V12 = 47rV (P+x(5i, ^2) + P-X(l - <52, 1 - S^)) . (26) 

The value of 6q^ associated to it is 



'^-'P-^7u-i))""'-Ul+»' (^^' 



These expressions can now be introduced into the general formula Eq. ([H]) to 
obtain the first solution 



M« = -^i^P^cj^dcui-62) - P-(r,duj{52: 

-t7cp[Ju\z + -)P^a^duj{z)+ J u\z+-)P+a^duj{z))+h.c.. (28) 

<5i 1-52 

The quantities u and u are the ones associated to the Q = 1 caloron. The 
analytic expressions needed to do the calculation were explicitly given in our 
previous paper (TTj. 



Now we investigate the other solution, associated to 5q^ = 0. One has to solve 
the homogeneous equation (IMI) for vanishing right hand side. A particular 
solution is given by 



iJuiz) = exp{-r3 / dz' AA{z')} = 0,(z) - t^Mz) ■ (29) 



Since AA{z') is constant at intervals, the integral in the exponent is trivial 
to perform. We leave the explicit form of 0s (-z) and 0a (-2) to the reader. It is 
interesting to point out nonetheless, that 0s (-2) is periodic in z with period | 
and 0a (2) antiperiodic. 

From the previous expression we can compute the corresponding self-dual 
deformation using Eq. ([8]). The result is given by 



542) = ^Jdz {u^{z + ^)(0s(^) + T,U^))a,duj{z)) + h.c. . (30) 



Again, the integration over z can be performed analytically using the formulas 
of our previous paper pjj . 

We have arrived to the general solution our problem. The adjoint zero-modes 
of the (self-dual) caloron which are antiperiodic in time are 

vl> = ^M«7m(I + 75)V^i + bA^;^h,iI + 75)V^2 , (31) 



where Va are arbitrary constant spinors and 5A^^' are given in Eqs.( !28l) -( !30l) . 

It is interesting to mention that the general investigation of the possible values 
of 6q^ has led to another solution having a fairly simple form. The expression 
of the left-handed Weyl spinor, \1>(^) = \l/^^Va, is: 

¥^^ = a.d.T'^aaV , (32) 



where a labels a colour component, cXa acts on the spin indices and V denotes 
an arbitrary constant 2-spinor. The functions Ta depend on the colour index 
as: T^ = T^ = —1/F and T^ = P+x + P-X- The function x is essentially 
the function with the same name given in Ref. [2|5] . Curiously this solution 
interpolates between the non-supersymmetric periodic adjoint zero-mode for 
mi = {Si = 0) and one of our antiperiodic solutions (Eq. fl28|) ) for mi = 1712 
{61 = 1/4). Using the formulas of the next section it can be proven that the 
solution is neither periodic non antiperiodic for other values of the mass mi. 



3 Properties of the solutions 



In this section we will investigate the general properties of the solutions found 
in the previous section. 



3.1 Periodicity in time 



Here we will explicitly verify the required antiperiodicity in time of our general 
solution. In our gauge the caloron vector potential satisfies 

4°)(xo + 1) = e*"¥^4°)(xo)e-^^ (33) 



Thus, the required antiperiodicity of the adjoint zero-modes amounts to: 

. m.1 To . JTil TQ 

^(xo + 1) = -e'^^^{xo)e-'^^ (34) 

This property follows easily from the form of our solutions and the periodicity 
behaviour of u: 

u{z, xo + 1) = e*2""M(z, xo)e-'^^ (35) 

and an identical relation for uj and q. 

3.2 Far-field limit and Normalization 



The reader might question whether our general solution Eq. fl3T|) is nor- 
malizable. One can investigate the behaviour at points whose distance to 
the location of the constituent monopoles (ri and r2) is much larger that 
(3 and that vrp^. For the unequal mass case the zero-mode density goes to 
zero exponentially as e"'^™^"™'^)''^. The equal mass case {mi = m2 = vr) 
is more subtle since both solutions decay in power-like fashion. The non- 
homogeneous solution Eq. (!28l) coincides with the additional solution Eq. (J32l) 
in this case. In the limit under consideration x goes to zero exponentially and 
F = {ri + r2 + 7rp^)/(ri + r2 — vrp^). Thus the density behaves as 1/r^. 

An alternative approach to normalizability of the solutions is to compute the 
norm of the solutions. In fact there exist a general formula [2][T9] which allows 
one to compute the norm and the scalar products of the solutions in terms of 
Nahm-data directly. This is also useful in checking if the real dimensionality of 



10 



the space of solutions is 8 (4 complex dimensions, 2 quaternionic dimensions), 
as indicated by the index theorem. Using this formula we obtain 



5A^^ 



p- 



2 



Air' + SvrV^ {62 - 61) , (36) 

p ^ sinh(4.V^0 ^ (,^ _ ,^) ,,,h(4.V^^i) , (37) 

(M«, Mf ) = 2n'p e-2-V^^i (5^ _ 5^) . (38) 

3.3 Profile of the zero-mode density 



In this subsection we will describe the qualitative properties of the zero-mode 
densities. For that purpose we developed two independent programs to draw 
these profiles. Both programs give matching results. In Fig. [1] we give the 
contour plot in a z-y plane of the solution 5A^^^ (top) and an orthogonal 
CP-pair (bottom) for p = 1 (giving an intermediate size caloron separation) 
and two representative values of the masses. The z axis is the line joining 
the constituent monopoles and is represented horizontally. The vertical axis 
denotes the y axis (the density is axially symmetric). 

For the equal mass case [mi = 7722 = tt) the mode following from Eq. (|30l) has 
an approximately constant higher density along the line joining both calorons 
(top left). This can be interpreted as a string. In contrast, the other solution 
associated to Eq. fl28l) has a region of small density located along the line 
joining the two monopoles (bottom left). As the masses become unequal, the 
most massive monopole dominates the densities. The right contour plots show 
the situation for 61 = 0.23. 



3.4 Limiting cases 



The caloron is an interesting solution which interpolates between the gauge 
potential of an instanton and that of a BPS monopole. It is interesting then 
to see how our antiperiodic zero-modes behave in these extreme cases. We 
will first concentrate in the situation corresponding to the trivial holonomy, 
Harrington-Shepard, caloron: 5i = 0. In that case one of the constituent 
monopoles is massless and pushed to infinity. The p parameter of the so- 
lution does no longer control the separation between the monopoles but is 
still a free parameter. For small p the HS caloron approaches an ordinary, 
zero temperature, instanton. From our general formulas, it is easy to check 
that in that limit and close to the center of the caloron the time periodicity 
becomes irrelevant and the two zero mode CP-pairs approach the periodic 
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Fig. 1. Contour plots of the density of the two antiperiodic zero modes in the y — z 
plane. Constituent monopoles are localized at y = and separated along the z axis 
which is drawn horizontally. Left: For mi = m2 = vr and /J = 1. Right: For 5i = 0.23 
and p = 1. 

zero modes of the instanton. In the opposite, p ^ oo, limit the HS caloron 
becomes a BPS monopole with time independent action density. Despite the 
time independence of the background there are still 4 non-trivial antiperiodic 
zero-modes. They can be easily derived from Eqs. fl28|) and (!30|) by taking 
the appropriate (5i = and p — > cx) limits. Up to a gauge transformation we 
obtain: 



MW'(i) = %",. irp (e\(a:) E^^ir) - el(x) K(r] ] + h.c. 



MfW = ^ 



C-O \ iXj 



^^W+e1(x)4(r))+h.c. 



(39) 
(40) 



where SA^^f^ is directly derived from Eq. (|30l) and SA^^^' is the combination of 
Eqs. (128|) and (130!) orthogonal to 5A^^\ In the expression above, E^^ is the 
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electric field of the BPS inonopole: 

sr(.) = -,£!(^P„.-,MMp-. (41) 



and we have introduced the time independent quantity: 

2cosh(7rrj v ^^ f / 



with g{u) = M/sinh(M), g'{u) its derivative with respect to u, and P^ = 
(a^ ± n(j^j_h)/2, n = XiTi/r. The antiperiodicity of the solution is encoded in 
the time dependent quaternions ef and el defined through: 

g-iTTx^s^ = i{el{x) + ielix)) . (43) 



For non-trivial holonomy there is also an interesting limit in which the caloron 
solution tends to the BPS monopole. It corresponds to making the separation 
of the constituent monopoles tend to infinity (p -^ oo). Our adjoint zero- 
modes lead to those of the BPS monopole if the appropriate limit is taken 
(ri << Tcp^, p » 1). For example, for the equal mass case (mi = m2 = vr) 
the first solution Eq. fl25]) follows quite simply by applying the appropriate 
hmit to Eq. (152]) . Computing the density we obtain: 

2{h'^{7rr) + 1 - 2h'{TTr) cos0) + /(vrr) + g'^{7Tr) + 2g{7ir)g'{7Tr) cos (44) 



where h'{u) is the derivative of h(u) = Mcoth(M). This profile has axial sym- 
metry depending explicitly on the azimuthal angle 6. Notice also that the 
solution is non-normalizable. 



3.5 Comparison with numerical results 



We have crosschecked our results with a direct evaluation of adjoint zero- 
modes on the torus obtained by lattice methods using Neuberger's overlap 
operator [17J in the adjoint representation. One expects that the spatial pro- 
file of the torus solutions approaches our analytical formulas as the box size 
becomes much larger than all scales of the problem ( (3 and vrp^). To make 
a quantitative comparison we computed the zero-mode density integrated in 
time along the line x = y = joining both constituent monopoles. It is not 
possible a priori to construct numerical zero-modes with a prescribed value of 
p and ^1, although some tuning is possible |18j. For the numerical comparison 
displayed in Fig. [2] we slightly tuned by hand these parameters to improve the 
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Fig. 2. Comparison between numerical (circles) and analytic (lines) zero modes for 
p = 0.79 and 5i = 0.172. We display the density of the zero modes, integrated in 
time, along the line joining the two monopoles. 

agreement (p = 0.79, 5i = 0.172). A technical point which one has to address 
is how to guarantee that the same linear combinations are selected for the 
numerical and analytical data. We chose to define the two linearly indepen- 
dent modes by imposing that at the center of mass {x = y = z = 0) one has 
maximal and the other minimal density (integrated over time). 



4 Conclusions 



In this paper we have obtained analytic formulas for the zero-modes of the 
Dirac equation for gluinos in the background field of Q=l SU(2) calorons with 
antiperiodic boundary conditions in thermal-time. Our formulas are valid for 
non-trivial holonomy as well as for the Harrington-Shepard caloron and include 
as a limiting case those of BPS monopoles. The solutions have finite norm and 
decay exponentially with distance if the masses of the constituent monopoles 
differ. Their density profile contrasts with the case of periodic zero-modes. For 
example, as the monopoles are pulled apart the density does not decouple into 
independent lumps centered at the monopoles, but rather describes a string 
joining the monopoles. Nonetheless, the number of normalizable zero-modes 
matches in both cases. 

Our work has methodological interest since our approach is applicable to other 
cases including the extension to SU(N), and might be instrumental in finding 
formulas for calorons of higher charge. From a physical viewpoint our work 
provides a first step towards a semi classical study of N=l SUSY Yang-Mills 
at finite temperature. There are interesting issues at stake such as that of 
supersymmetry breaking at finite temperature, which has been a subject of 
debate since early times p0]-[26]. It is our intention to address these questions 
in future work. 
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